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ABSTRACT 


The  design  and  evaluation  of  optimum,  time -invariant  filters 
(the  classical  approac  i  to  prediction  problems)  are  subject  to  several  limi¬ 
tations.  Because  of  these  limitations,  an  adaptive  system  was  considered. 

This  report  describes  a  filtering  system  which  can  adjust  to  changes  in  either 
signal  or  noise,  thereby  overcoming  many  difficulties  of  time -invariant 
filtering.  The  small  amount  of  required  computational  time  to  update  the  filter 
weights  is  another  important  feature  of  the  scheme  presented. 

In  this  report,  an  adaption  algorithm  applied  to  a  simple  time- 
series  model  was  studied.  For  the  case  of  stationary  data,  a  tradeoff  between 
the  adaption  rate  and  the  mean -square -error  performance  of  the  filter  exists. 
The  adaption  rate  is  inversely  proportional  to  the  convergence  parameter  X, 

while  the  mean-square -error  is  directly  proportional  to  X.  For  the  case  of 
nonstationary  data,  a  tradeoff  between  adapting  too  slowly  and  adapting  too 
rapidly  exists.  The  optimum  rate  of  adaption  appears  to  be  approximately  10 
times  faster  than  the  average  time  rate  of  change  in  the  input  data  statistics. 
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SECTION  I 

THE  ADAPTIVE  ALGORITHM 


The  classical  approach  to  prediction  problems  is  to  design  and 
evaluate  optimum,  time-invariant  filters.  Unfortunately,  these  filters  are 
subject  to  several  limitations  which  restrict  their  usefulness.  Numerical 
determination  of  the  optimum  filter  is  often  involved  and  costly.  Even  after 
the  optimum  filter  is  determined,  it  *s  optimum  only  for  the  signal  and  noise 
models  from  which  it  was  designed.  A  change  in  the  signal  or  noise  model 
requires  a  new  derivation,  which  frequently  is  considerably  difficult  and 
expensive.  This  report  describes  a  filtering  system  which  can  adjust  to 
changes  in  either  signal  or  noise  and  thus  overcome  many  difficulties  of  time- 
invariant  filtering.  The  small  amount  of  required  computational  time  to  up¬ 
date  the  filter  weights  is  another  important  feature  of  the  scheme  presented. 

In  general,  adaptive  systems  can  be  considered  the  combination 
of  two  systems  —  an  adjustable  "worker,  "  with  internal  parameters  which 
can  be  varied  in  value  upon  command  from  a  "supervisor. "  The  purpose  of  the 
supervisor  is  to  measure  the  filter's  performance  and  to  output  the  necessary 
adjustment  to  the  worker  (Figure  1).  Of  particular  interest  are  adaptive 
devices  which  perform  a  linear  operation  on  the  input  data.  The  adaptive 
process  can  be  viewed  then  as  the  search  for  a  weight  vector  that  "optimizes" 
the  response  of  the  filter. 

There  are  many  ways  to  measure  how  well  the  system  performs 
its  intended  task.  One  of  the  most  common  (the  one  which  shall  be  adopted  in 
this  paper)  is  the  criterion  of  least-mean- square-error,  the  criterion  used  in 
the  Wiener  theory  of  signal  extraction  and  prediction.  1  Not  only  is  the  least- 
mean- square-error  criterion  convenient  for  comparison  to  classical  time- 
invariant  filtering,  but  it  also  lends  itself  quite  easily  to  an  adaptive  algorithm, 
which  shall  now  be  considered. 
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Figure  1.  An  Adaptive  Processor 


A.  DERIVATION  OF  ADAPTIVE  ALGORITHM 

The  adaptive  algorithm  can  be  derived  by  onsidering  the 

Wiener  theory  for  sampleu-data  systems.  The  theory  is  concerned  with 

finding  the  linear,  time-invariant,  sample-point  operator  H(t),  which  weighs 

the  input  time- series  jx  }  in  such  a  way  that  the  output  of  the  filter  at  time 
'c  •  1 

t,  S^,  is  the  best  mean-square  approximation  to  the  desired  output  S^. 

Assuming  that  the  impulse  response  of  the  filter  is  of  the  form 
M+L-l 

h(t)  "  2  aj  6  (t- 3  At)  (1) 

j=L 

the  mean- square-error  can  be  written  as 


C  = 


lim 

T  «•  eo 


<f  S 


t=l 


lim 
T  -*  eo 


M+L-l 


a.  x. 

J  t-j 


(2) 


The  filter  weights  aj  which  minimize  Equation  2  are  those 

weights  for  wluch  the  gradient  of  7?  equals  0.  Solving  for  the  a.,  the 

2  .  J 

following  matrix  equation  is  obtained: 

$ 

£•  a  =  iE  (3a) 


In  this  report,  a  double  bar  under  c.  literal  expression  represents  a 
square  matrix,  and  a  single  bar  represents  a  column  matrix. 
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For  the  single- channel  problem  of  predicting  L  points  ahead 
with  an  M-point  filter,  the  elements  R,  a,  _cp  can  be  written  as 


R  = 


r  (0)  r  (1)  . 

x  x 

r  (1)  r  (0)  . 

x  x 


•  •  • 


r  fM-1) 
x 


•  #  •  • 


r  (M-l)  .  .  .  .  r  (0) 

X  x 


(3fc) 


and 


r  > 

a  t  I 


a  =  < 


M+L- 1 

v.  -/ 


(3c) 


r  r  (L)  ^ 

x 


s.  =  \ 


r  (M+L- 1) 

L  x  J 


(3d) 


where  r  (j)  is  the  j**1  lag  correlation  value  for  the  data. 
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Rather  than  solving  Equation  3a  for  the  filter  a  directly  by 
inverting  the  R  matrix,  an  iterative  method  can  be  used  to  obtain  successive 
approximations  to  the  optimum  filter  weights.  In  the  technique  considered, 
a  new  filter  matrix  a  (k+1)  is  determined  from  the  previous  filter  matrix  a  (k) 
by  traveling  in  the  direction  which  maximizes  e  change  of  the  mean-square- 
error  for  a  given  change  of  the  filter  matrix,  ‘-’his  direction  is  given  by  the 
gradient  vector  of  the  mean-square-error: 


The  optimum  step  length  X  down  the  gradient  for  each  iteration 

can  be  found  in  a  systematic  way  from  iteration  to  iteration.  Fortunately  for 

computational  considerations,  \  can  also  be  held  constant  and  the  method 

3 

will  still  converge.  Thus,  our  scheme  can  be  written 

a  (k+1)  =  a  (k)  -  y  1  *  (4) 

a  =  a  (k) 

A  drawback  of  the  iterative  procedure  derived  in  Equation  4 
is  that  it  requires  knowledge  of  the  gradient  of  the  mean-sqaare-error  and 
this  is  seldom,  if  ever,  known  in  practice. 
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Several  authors 


4,5 


have  shown  that  the 


preceding  algorithm 


can  be  modified  by  using  the  instantaneous  estimate  of  V 
given  by 


a  =  a  (k) 


The  adaptive  algorithm  then  becomes 

a(k+l)  =  a(k)  +  \c  x  (5) 

This  scheme  can  be  generalized  to  an  NC-channel  problem  by 
treating  the  data  matrix  as  being  composed  of  data  from  all  NC  channels. 
The  filter  matrix  a(k)  represents  the  NC  filters  applied  to  the  data  vector 
Figure  2  shows  this  filter  system  schematically.  The  superscript 
represents  the  channel,  and  the  subscript  represents  the  time  point. 

B.  TIME-CONSTANT  CONSIDERATIONS 

It  is  desirable  to  be  able  to  relate  the  parameter  \  to  the 
amount  of  data  being  used  to  design  the  set  of  filter  coefficients  at  each 
point.  We  propose  to  do  this  in  the  following  manner. 


o* 
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Figure  2.  An  Adaptive  Multichannel  Filter  System 


We  define  the  sequence  \  a  '  }  by 


*  ^k  +  xk  £k-xk 


where 


"  T 


v  V 

-k  —  k 


and  is  the  adaptive  filter  matrix  at  time  k.  The  sequence  can  be 
considered  the  “false  gain"  set  of  coefficients,  since 


2k 


r-k  ‘  +  *k) 2  0 


We  consider  the  sequence 


~  -i/T  , 

2-.  e  a 


i=  1 


-k-i 


”k  v  -i/T 
L,  e 
i=  1 


/,  -1/T\  ,  -1/T  . 

are  jin  +  6  ^k-i 


*  in +  (*-e’i/T)Uk-i  -  in) 


T  denotes  transpose 
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If  we  let  k  then 


4+i  -  4  +  4  (*  *  e’  0 


-  1/t\ 

e  )  c  ,  x. 


where 


=  4  +  Vk4 


4  *  4'  (‘  -  ■'  ) 


If  X'  is  constant  in  time,  note  that  this  algorithm  is  the  adaptive  algorithm 

K 

derived  earlier,  with  a  time  constant  T. 

In  general,  X.^  varies  in  time;  the  degree  of  variation  affects 
the  time  constant  of  the  adaptive  algorithm  in  Equation  5;  however,  a 
"ball  park"  estimate  of  T  can  still  be  obtained.  Take  X.  '  to  be  given  by 


tk2kT] 


M  •  NC  •  r  (0) 
x 


where 


M  =  number  of  filter  points 


NC  =  number  of  channels 


r  (0)  =  zero-lag  correlation  of  the  data 
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The  ratio  between  X '  and  X  used  in  Equation  5,  then, 
will  give  an  estimate  of  t  ;  i.e. , 


h.  -  _ X _  -  1/t 

X.  "  M  •  NC  •  r  (0)  " 
k  x'  ' 


(6b) 


Appendix  D  shows  that  time  constant  t  can  be  associated  with 
an  "equivalent"  boxcar  section  of  data  of  length  N  by  the  relation 


N 


1  -f  e 
1  -  e 


1/t 

1/t 


(6c) 
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For  the  multichannel  adaptive  filter,  the  behavior  of  the 
algorithm  becomes  obscured  by  the  complexity  of  the  system.  For  this 
reason,  the  adaptive  process  is  examined  for  a  simple  single- channel 
prediction  problem.  The  results  of  this  investigation  give  some  intuitive 
feeling  for  the  behavior  of  this  scheme  for  more  complicated  filtering 
problems. 

A.  CASE  1,  STATIONARY  TIME  SERIES 

One  of  the  simplest  time  series  to  which  the  adaptive  filter 
could  be  applied  is  a  first-order  autoregressive  series.  The  series  |xt| 
is  defined  by  the  recursive  relation 


xt  =  axt-i  +  nt  <7> 

where  the  series  |n^|is  a  sequence  of  independent  Gaussian  random  variables 
of  zero  mean  and  unit  variance.  The  parameter  a,  which  must  be  less  than  1 
in  magnitude  for  stability,  determines  the  coherence  (predictability)  of  the 
time  series. 

For  the  case  where  a  is  a  constant  in  time,  the  process 
defined  by  Equation  7  is  stationary.  The  optimum  Wiener  prediction 
filter  is  given  by 

% 

h(t)  =  a  6  (t-  1) 

with  a  mean-square-error  of  unity.  Thus,  we  can  analyze  the  1-point 
adaptive  filter  by  comparing  the  filter  weight  to  the  optimum  time -invariant 
filter  weight  a  and  the  mean-square-error  of  this  process  to  unity. 

6  is  the  impulse  function. 
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By  starting  the  adaptive  filter  weight  a  at  0  and  letting  the 
filter  adjust,  we  can  measure  the  time  constant  of  adaption.  On  the  basis 
of  theoretical  considerations,  a  reasonable  estimate  of  the  time  for  initial 

transients  to  decrease  to  1/e  of  their  original  value  is 

* 

1 

\r  (0) 

x 

where  r^O)  is  the  zero-lag  correlation  of  the  data  and  X  is  the  convergence 
parameter.  The  behavior  of  a^  has  been  studied  experimentally  by  simulating 
the  adaptive  process  on  an  IBM  7044  computer  at  Texas  Instruments  Incor¬ 
porated.  Results  of  this  investigation  for  various  values  of  X  for  a  fixed  a 
are  shown  in  Figure  3  where  the  dashed  line  is  the  theoretical  curve  for  the 
expected  value  of  a We  see  that  the  theoretical  and  experimental  curves 
are  in  close  agreement 

The  area  of  primary  interest  is  the  L.MS  adaptive  filter's 
steady-state  behavior.  To  compare  the  steady-state  performance  of  the 
adaptive  system  with  that  of  the  classical  time -invariant  filter,  it  is  con¬ 
venient  to  introduce  a  "misadju**ment"  parameter  M.  The  misadjustment 
M  is  a  dimensionless  quantity  defined  by  the  ratio 


MSE  -  MSE 

a  o 

M  '  MSE 

o 

(8) 

where 

MSE  * 

a 

mean-square-error  of  output 

CL 

adaption  filter 

MSE  = 

f) 

mean- square- error  of  output 

classical  time -invariant  filter 

* 

See  Appendix  A  for  the  theoretical  derivation. 
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The  misadjustment  coefficient  M  has  been  experimentally 
determined  for  the  autoregressive  time  series  defined  by  Equation  7. 

Results  for  three  values  of  a  as  a  function  of  \  are  plotted  in  Figure  4  where 
the  dotted  lines  show  the  theoretical  curves  for  the  misadjustment  coeffic¬ 
ient  as  a  function  of  The  filter  weight  is  started  at  the  proper  value 
a^  =  a,  and  the  mean-square-error  is  computed  over  1000  points.  The 
curve  for  a  =  0.2  (Figure  4)  falls  below  the  0-percent  line  because  of  the 
statistical  fluctuations  in  the  data.  Had  the  mean- square-error  been  com¬ 
puted  over  an  infinite  length  of  data,  this  curve  would  have  been  slightly 
above  the  0-percent  line. 

B.  CASE  2,  NONSTATIONARY  TIME  SERIES 

Of  practical  importance  is  the  behavior  of  the  adaptive  filter 
system  when  the  input  time  series  is  nonstationary.  Data  used  for  this 
experiment  are  first-order  time-varying  autoregressive  Beries 

*,  =  at*t-i  +  nt 

where  n  is  as  previously  defined  and  a  =  A  +  B  sin2TTft.  A  is  the  mean 
t  t 

value  of  a^,  B  is  the  amplitude  of  the  sinusoidal  deviation  about  the  mean, 
and  f  is  the  frequency  of  the  deviation. 

The  experimental  results  of  this  investigation  are  shown  in 
Figure  5.  Figure  5a  plots  the  mean-square-error  of  the  adapti\e  filter  as 
a  function  of  \  for  three  frequencies.  The  autoregressive  data  are  generated 
with  at  =0.4  sin  2TTft.  The  mean-square-error  is  calculated  over  5000 
points.  Figure  5b  plots  the  same  set  of  curves  for  a  =  0.  5  +  0.4  sin2nft. 


See  Appendix  B  for  the  theoretical  derivation. 
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Figure  5.  Mean -Square -Error  as  a  Function  of  X  for  Three  Frequencies 
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As  might  be  intuitively  expected,  the  optimum  value  of  X.  ir  creases  with 
the  increasing  time  variability  of  the  data  (increasing  frequency  of  0^). 
Also,  on  the  basis  of  the  results  with  the  stationary  data,  one  would  expect 
the  time-varying  data  with  the  largest  mean  value  of  a  (largest  average 
power)  to  be  the  most  sensitive  to  changes  in  X..  This  is  indicated  by  the 
results. 


The  performance  of  the  adaptive  filter  appears  to  be  optimized 
when  the  filter  can  adjust  to  the  data  at  a  rate  10  times  faster  than  the  rate  of 
change  in  the  data  statistics.  In  other  words,  the  optimum  value  of  X.  can  be 
found  from  the  empirical  relation 


10 


X  r  (0) 
opt  x  avg 


(9) 


where 


da 

1 

~dt 


avg 


=  4Bf  is  the  average  rate  of  change  of  the  autoregressive 


coefficient,*  r  (0)  is  the  average  power  in  the  trace,  and  X.  is  the 
x  avg  6  r  opt 

optimum  value  of  the  convergence  parameter.  In  Table  1,  the  X  given 

opt 

by  Equation  9  is  compared  with  the  true  X  ^  for  each  of  the  six  experiments. 

Without  additional  theoretical  results,  further  experimental 
study  of  the  adaptive  algorithm  as  applied  to  nonstationary  autoregressive 
time  series  is  hardly  merited.  It  would  be  of  interest  to  study  its  behavior, 
however,  if  the  convergence  parameter  were  allowed  to  change  in  some 
optimum  manner  as  the  filter  "tracked"  the  data.  Comparison  of  the  results 
of  such  an  investigation  with  those  presented  in  this  paper  would  give  some 
indication  as  to  whether  this  approach  would  significantly  improve  the  per¬ 
formance  of  the  adaptive  filter. 


See  Appendix  C  for  the  derivation. 
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AND 


SECTION  III 
CONCLUSIONS 


i  or  this  report,  we  have  studied  an  adaption  algorithm  applied 
to  a  simple  time-series  model. 

For  the  case  of  stationary  data,  it  is  shown  that  a  tradeoff 
between  the  adaption  rate  and  the  mean-square-error  performance  of  the 
filter  exists  and  that  the  adaption  rate  is  inversely  proportional  to  the  con¬ 
vergence  parameter  while  the  mean- square- error  is  directly  proportional 
to  \. 

.For  the  case  of  nonstationary  data,  a  tradeoff  between  adapting 
too  slowly  and  adapting  too  rapidly  exists.  The  optimum  rate  of  adaption  ap¬ 
pears  to  be  approximately  10  times  faster  than  the  average  time  rate  of  change 
in  the  input  data  statistics. 
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APPENDIX  A 


THEORETICAL  BEHAVIOR  FOR  EXPECTED  VALUE  OF 
ADAPTIVE  FILTER  WEIGHT 

For  a  1 -point  prediction  filter,  the  adaptive  algorithm  of 
Equation  5  becomes 

a  (k  4  i)  ^  a(k)  +  X‘kxk  l  (A-l) 

where 


€k  =  Xk  "  a(k)Vl 


(A-2) 


Substituting  Equation  A-2  into  Equation  A-l  and  rearranging 


terms,  we  obtain 


a  (k  +  1) 


’  (‘  '  Xxl>) 


a(k)  +  X  xk  xkl 


(A-3) 


The  relation  between  the  k  +  1  iteration  and  the  initial  filter 
weight  a  (1)  is  obtained  by  recursive  substitution  into  Equation  A-3,  resulting 


a (k  +  1) 


=  (1-XxLi)a<k)  +  Xxkxk-i 
*  (’  ■  Xxk-i)  (*  -  a(k  -  11  +  xVi  Vz] 
+  Xxk  xk-i  *  (‘  *  Xxt-i)(‘ '  Xx k-2) a ^ 

+  Xxk  xk-l  +  X(‘  "  Xx  k-l)  xk-l  xk-2 


+  Xxkxk 


(A-4) 


“[A  ('  ■  Xxi-i)]a<i)  +  xi5I  [xixi-i  j?+1  (■  •  Xxj-i) 
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using  the  convention 


1  for  p  >  k. 


Consider  Equation  A-4  for  the  case  where 
compared  with  one  for  all  j.  We  can  then  approximate  a 


Xx  .  ,  is  small 
J-1 

(k  +  1)  by 


,  2 

“ 

2 

a  (k  +  1)  ss 

k  -  Xx 

n  1-1 

n  e 

k 

a(l)  +  X  £ 

X,  x.  , 
l  l-l 

k  -Xx  .  , 

n  e  J‘‘ 

i=  1 

i=  1 

Ci. 

II 

»-*• 
+ 
H- » 

. 

_ 

J 

(A-5) 


k  -Xx^ 

n  e  l_1 
i=l 

k 

a  (1)  +  Z 

k  -X*2  , 

n  „  J’1 

i= 1 

V  / 

+ 

•H 

II 

•n 

* 

- 

J 

The  expected  value  of  a  (k  +  1)  is  approximated  by 


(A-b) 


r  i 

’k  -Xr  (0)' 

k 

7  -Xar  (0)\  k  -Xr  (0)" 

a  (k+  1)  « 

n  e  x 

a  (1)  +  E 

(l  -  e  )  n  e  X 

L  J 

i=l 

i=l 

\  /j=i+l 

assuming  that  the  data  are  stationary  where  r^(0)  is  the  zero-lag  correlation 
of  the  data  and  ar^(O)  is  the  first-lag  correlation  value. 

Simplifying  Equation  A-6,  we  obtain 


E  |a  (k  +  1)J  «a(i) 


-Xk  r  (0) 
x 


Xar  (0) 


AV**  1VI  r  *1 

1  .  e  x  f  "Xk  r  (0)1 

+  T~^r  t-  J 


-Xk  r  (0) 

X 

»  a  (1)  e  +  a  1 1  -  e 


r  -Xk  r  (0)  “I 
x  .  x 

a  1  -  e 
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This  expression  shows  that,  on  the  average,  the  initial  adap¬ 
tion  transient  should  decay  with  a  time  constant  of  - — per  iteration. 

\r (0)  r 

Also  there  is  the  result  that,  in  the  limit  as  k  approaches  infinity,  the  ex¬ 
pected  value  of  our  filter  weight  is  the  optimum  value  a;  i.  e. , 


lim  E  [a  (k  +  1)1  =  a 
k  «•  <=>  L  J 
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APPENDIX  B 


DERIVATION  OF  MISADJUSTMENT  COEFFICIENT  FOR 
1 -POINT  ADAPTIVE  FILTER 

In  the  text,  we  define  the  mis  adjustment  of  an  adaptive  filter 


MSE  -  MSE 
_ a _ o_ 

MSE 


where 


MSE  =  mean -square -error  for  the  output  of  the 
adaptive  filter 

MSE  =  mean -square -err or  for  the  output  of  the 
classical  time -invariant  filter 

We  now  evaluate  thic  quantity  on  the  basis  of  theoretical 
considerations  for  the  1 -point  prediction  *ilter. 


given  by 


The  error  term  of  the  adaptive  filter  for  the  kth  output  is 


e,  =  x  "  a  (k)  x. 


=  xk-(ixk-i- 

where  is  the  difference  between  the  adaptive  filter  weight  and  the 
optimum  time -invariant  filter  weight.  The  expected  value  for  the 
mean -square -error  for  the  k**1  output  is 


(B-l) 


E  E  [(xk_  aitk-i)2-24kxk-i(*k-'oik-i)+  Ak  *k-i] 

'  -  MSEo  +  *  [»* 


(B-2) 


B-l 
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To  evaluate  the  second  term  in  Equation  B-2,  we  can  obtain,  by  recursive  substitution 


B-2 
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APPENDIX  C 


DERIVATION  OF  AVERAGE  RATE  OF  CHANGE  FOR 
TIME-VARYING  AUTOREGRESSIVE  COEFFICIENT  afc 

In  the  text,  we  define  the  time -varying  autoregressive  coeffi¬ 
cient  to  be  of  the  form 


a  =  A  +  B  sin  2irft 


(C-l) 


where  A  is  the  mean  value  of  afc,  B  is  the  amplitude  of  the  sinusoid^1  devia¬ 
tion  about  the  mean,  and  f  is  the  frequency  of  the  deviation.  The  average 
rate  of  change  of  is  given 


da 


dt 


=  2rrfB  |sin  2nft] 
avg  avg 


(C-2) 


where  |  |  indicates  the  average  absolute  value.  Clearly, 

avg 


sin  2nft 


tt/2 

sin  xdx  =  — - 
avg  TT  J  TT 

o 


-u 


(C-3) 


Substituting  the  above  value  into  Equation  C-2,  we  obtain  the 
result  stated  in  the  text: 


dt 


=  4Bf 


(C-4) 


avg 
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APPENDIX  D 


DERIVATION  OF  "EQUIVALENT"  BOXCAR  LENGTH  N  ' 

WITH  EXPONENTIAL  TIME  CONSTANT  t 

In  the  text,  v  e  define  the  exponential  weighting  scheme  as 


ao  .  / 

^  “1/T 


E  W  T  > 

s  a  . 

i=i _ (.  -i/T  , 

Ik  "  F  "T7t  “  V  ‘  6  /  E  e  -k-i 

itl e  i=1 

The  variance  of  the  random  variable  £  is  given  by 

[«,]  ■  «  [is,  ■  « [t,l )‘] 

■  (■  *[,?,  £  ■  ■[=;., 

If  the  a^  ^  are  uncorrelated  time- stationary  random  variables. 


E  fev  -  E 


W)1 


■(■  I,  *“'T  -u  ] 


1  -  e 


■1/t 


1  +  e-1/T 


W 
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It  is  a  well-known  result  that  the  variance  of  the  average  of 

•  •  w  6 

N  random  variables,  each  with  the  same  variance,  is  given  by 


aN  = 


1  2 
-rr  a 
N  £ 


Making  the  identification 

-1/t 


1  _  1  -  e 

N  =  .  . 

1  +  e 


-T7t 


we  have  the  result  that  the  time  constant  Tmay  be  associated  with  an 
"equivalent"  boxcar  section  of  data  of  length 

i  “1/t 

1  +  e 

N  - - rrr 


1  -  e 
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